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Foreword 

The  necessity  of  determining  the  liquid  water  content  along  a  trajectory  or 
iine-of-sight  recurs  ever  more  frequently  as  more  stringent  performance  specifica¬ 
tions  are  placed  on  navigation,  communication  and  weapon  systems  operating  at 
ever  smaller  wavelengths  in  the  microwave  spectral  region.  Scientists  at  AFGL 
(and  its  predecessor,  AFCRL)  have  long  had  an  interest  in  the  effects  of  clouds 
and  precipitation  on  microwave  propagation,  and  have  done  much  theoretical  work 
relating  observed  physical  characteristics  of  hydrometeors  to  attenuation  and  back- 
scatter.  1,2  However,  it  was  not  until  1971  that  the  opportunity  arose  to  combine 
aircraft  measurements  of  cloud  and  precipitation  with  ground-based  radar  observa¬ 
tions  on  a  systematic  basis.  At  the  request  of  SAMSO  (now  BMO),  AFGL  scientists 
provided  environmental  data  for  ballistic  missile  test  firings,  first  at  NASA's 

Wallops  Flight  Test  Center  (WFTC)  in  Wallops  Island,  Virginia,  and  later  (1973) 
w  3 

at  the  Army's  Kwajalein  Missile  Range  (KMR)  in  the  Marshall  Islands.  The  task 

as  originally  stated  was  to  determine  the  liquid  water  content  of  water  and  ice 
hydrometeors  encountered  by  a  vehicle  passing  through  the  atmosphere  at  hyper¬ 
sonic  speeds.  This  was  later  extended  to  include  the  determination  of  particle 
spectra  of  all  hydrometeors,  both  solid  and  liquid. 

The  requirements  of  the  weather  erosion  program  were  more  rigorous  than  any 
imposed  by  the  microwave  studies.  While  the  latter  could  be  met  by  an  integrated, 
or  average,  value  of  the  critical  parameter,  the  application  to  nosecone  erosion 
required  point  by  point  values  along  the  trajectory.  This  is  because  the  erosion  is 
a  function  of  missile  velocity  and  air  density,  or  altitude.  It  is  not  possible  to 
make  aircraft  measurements  along  the  expected  missile  trajectory  or  at  selected 
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altitudes  close  to  the  time  of  reentry.  Therefore,  it  was  determined  that  the  best 
way  to  meet  this  requirement  was  by  the  radar  measurement  of  the  radar  reflectivity 
factor  along  the  reentry  trajectory.  The  problem  then  became  the  determination 
of  accurate  and  precise  relations  between  the  observed  reflectivity  and  the  point 
values  of  liquid  water  content.  These  relations  were  established  by  combined  radar 
and  aircraft  measurements  taken  before  and  after  the  missile  reentry. 

Cloud  and  precipitation  particle  size  spectra  measurements  taken  aboard  an 
aircraft  are  characterized  by  small  sampling  volumes  and  loss  of  information  at 
larger  sizes  because  of  the  reduced  probability  of  detecting  these  fewer  particles 
and  also  because  of  instrument  truncation.  Radar  measurements,  on  the  other  hand, 
give  only  averages  of  the  summation  of  the  particle  backscatter  over  a  large  volume, 
and  do  not  record  contributions  from  the  lower  end  of  the  size  spectrum  because  of 
sensitivity  limitations.  Both  sampling  volume  and  sensitivity  are  functions  of  radar 
range.  The  problem  inherent  in  relating  aircraft  data  and  radar  measurements 
consumed  much  of  the  effort  devoted  to  this  task.  Some  of  the  solutions  arrived  at 


are  described  in  the  first  two  sections  of  this  report.  Section  1  describes  this 
method  of  extrapolating  an  instrumentally-trunc&ted,  aircraft -obtained  distribution. 
The  introduction  of  the  parameter  k  and  its  use  in  converting  from  aircraft  measure¬ 
ments  of  the  hydrometeors  to  radar  derived  values  of  M  is  contained  in  Section  2. 

The  scope  of  this  report  is  limited  to  methods  of  determining  M  which  includes 
the  mass  of  all  cloud  and  precipitation  particles,  whether  ice  or  water.  Other 
parameters,  saeh  as  average  particle  diameter,  while  important,  do  not  readily 
lend  themselves  to  radar  observation  without  previous  assumptions  concerning  the 
shape  of  the  size  spectrum. 

No  attempt  has  been  made  in  this  report  to  derive  and  evaluate  water  content  - 
radar  reflectivity  relations  within  the  melting  layer.  At  this  time  our  ignorance  of 
the  exact  nhysical  processes  and  their  interactions  taking  place  within  the  melting 
layer,  precludes  the  development  of  any  standard  procedure  for  remote  measure¬ 
ments  at  these  altitudes. 

Rosemary  M.  Dyer 
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TECHNIQUES  FOR  MEASURING  LIQUID  WATER 
CONTENT  ALONG  A  TRAJECTORY 

1.  Estimated  Distributions  From 
Instrumental ly  Truncated  Data 

Robert  O.  Berthe! 


1.1  INTRODUCTION 

The  specialized  instrumentation  used  aboard  the  AFGL  cloud  physics  aircraft 

includes  two  types  of  optic  tl  array  probes  manufactured  by  particle  Measuring 

Systems,  (PMS)  Inc. ,  of  Boulder.  Colorado.  Both  are  specifically  designed  to 

characterize  and  count  hyd'-ometeors.  The  one-dimensional  system  (PMS  1-D) 

classifies  hydrometeors  as  to  their  physical  size  in  one  dimension  and  also  measures 
4  * 

the  number  density.  The  t  ^-dimensional  system  (PMS  2-D)  provides  an  electron¬ 
ically  produced  shadowgraph  from  which  a  two-dimensional  area  can  be  determined 
and  particle  shape  inferred.  v 

Using  the  data  from  the  t-D,  the  liquid  water  content  (M)  of  the  sampled  hydro¬ 
meteors  can  be  determined  b-.- 


m  =  cL  N.  D? 


(Received  for  publication  11  March  1981) 

4.  Knollenlerg,  R.G.  (1870)  1  he  optical  array:  An  alternative  to  Scattering  or 

extinction  for  airborne  pa  rticle  size  determination,  J,  Appl,  Meteor. 
jh86-103. 

5.  Knollenbeig,  R.G.  (1976)  Three  new  Instruments  tor  ddud  physics  measure¬ 

ments:  The  2-D  spectrometer,  the  forward  scattering  spectrometer  probe 
and  the  active  scattering  aerosol  spectrometer.  Preprints  Internationa 
Conference  on  Cloud  Physics,  Amer.  Meteor.  Soc. ,  July  26-30,  1976. 
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An  investigation  of  the  instrument  imsratias  problem  has  resulted  in  a  method 
of  extrapolating  the  1-D  truncated  distributions  to  give  reasonable  approximations 
of  the  non-truncated  spectrum.  The  details  of  tMs  study  have  been  reported  else¬ 
where'  and  only  a  brief  description  of  the  technique  is  described  here. 

1.2  EXTRAPOLATION  OF  A  TRUNCATED  WSTROWnOlt 

The  problems  of  extrapolating  a  distribution  are  twofold.  Woich  method  is  to 
be  used  to  extend  the  distribution  beyond  the  upper  truncation  limit  and  what  is  the 
limit  of  the  extension,  or  the  maximum  particle  sire? 

When  addressing  the  first  problem,  it  is  logical  to  assume  that  the  form  of  the 
numbe r-deiisitv  distribution  will  remain  relatively  constant  beyond  the  physical  size 
limitation.  Therefore,  the  trend  of  the  distribution  should  be  established  from  the 
measured  portion  of  the  spectra.  It  was  found  using  the  1-D  large- snow  data  that 
in  the  majority  of  cases  the  distribution  is  adequately  described  by  an  exponential 


ot  tne  form 


»D  =  e 


Number  m  mm 


where  >'  is  the  intercept  of  the  distribution  function,  with  units  of  number  per  cubic 
meter  per  millimeter  band  with.  A  is  the  slope  of  the  distribution,  with  number  per 


size  CD  )  has  previously  oeen 
m 

distribution  by  an  expression  ol  tnc 


C  =  D. 


where  C  is  a  dimensionless  eons 
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These  findings  are  illustrated  in  Figures  1, 1  and  1. 2,  which  show  the  plotted 
number  density  vs  the  equivalent  melted  diameter  for  two  MC-130E  passes  taken 
on  4  July  1978  over  the  Kwajalein  Atoll.  The  hydrometeors  of  these  distributions 
are  classified  as  large  snow  and  have  maximum  physical  sizes  larger  than  the  1-D 
instrument  is  able  to  measure.  These  distributions  are  typical  of  the  experimental 
data  collected  by  AFQL  in  past  years,  exhibiting  a  systematic  decrease  in  the  numr 
ber  density  over  the  mid-  to  large-size  diameter  range.  They  also  demonstrate 
that  a  straight  line,  on  a  log-linear  plot,  can  adequately  describe  the  population  for 
the  major  part  (approximately  92%)  of  the  total  measured  size  range.  Deviation 
from  the  straight  line  occurs  at  or  about  the  junction  of  the  cloud  and  precipitation 
probes  with  the  cloud  probe  spectra  having  much  steeper  slopes.  The  contributions 
to  M  and  Z  from  the  precipitation  probe  data  (Mp,  Zp)  constitute  a  high  percentage 
of  the  measured  value  and  indicate,  at  least  in  truncated  situations,  that  the  M^and 
derived  from  the  cloud  probe  spectra  are  minor  contributions.  This  is  shown 
in  Table  1. 1. 


Figure  1. 1.  Number-Density 
Plot  of  Large  Snow  Lydrometeors 
as  Recorded  by  the  PMS  1-D  on 
the  MC-130E,  Pass  3  at  an 
Altitude  of  5. 1  km  Over  the 
Kwajalein  Atoll  on  4  July  1978 
(142  sec  average) 


When  the  data  from  the  precipitation  probe  of  the  1-D  instrument  are  considered 

by  themselves,  it  is  apparent  that  the  size  distribution  can  be  reasonably  described 

by  a  function  of  exponential  type  which  agrees  with  the  observations  made  by  other 
12  13 

investigators  for  snow.  ’  Although  this  occurrence  is  evident  in  the  majority 
of  cases,  occasional  sampling  levels  may  exhibit  number  density  distributions  that 
do  not  conform  to  an  exponential  fit  and  for  which  the  following  relations  do  not  hold. 

12.  Gunn,  K.  L.  b. ,  and  Marshall,  J.  S.  (1958)  The  distribution  with  size  of  aggre¬ 

gate  snowflakes,  J.  Meteorol.  15:452-479. 

13,  Sekhon,  R.  S. ,  and  Srivastava,  R.  C,  (1970)  Snow  size  spectra  and  radar 

reflectivity,  J.  Atmos.  Sci.  J2£:299-307. 
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Figure  1.  2.  Number-Density  Plot  of  Large 
Snow  Hydrometeors  as  Recorded  by  the  PMS 
1-D  on  the  MC-130E,  Pass  5  at  an  Altitude 
of  6.  6  km  Over  the  Kwajalein  Atoll  on 
4  July  1978  (212  sec  average) 


Table  1. 1  Comparison  of  the  M  and  Z  Contributions  Derived  From  the 
Measurements  of  the  Cloud  and  Precipitation  Probes  of  the  1-D 
Instrument 


Cloud 

Probe 

Precip 

Probe 

Total 

Percent 
Cloud  Precip 

..  -3 

M  gm 

Pass  3 

0.  0283 

0.3279 

0.3504 

6.4 

93.6 

Pass  5 

0.0244 

0.  3062 

0.  3306 

7.4 

92.6 

6  3 

Z  mm  m- 

Pass  3 

0.  054 

544.  1 

544. 154 

0.01 

99.  99 

Pass  5 

0.  040 

336.4 

336.440 

0.01 

99.  99 

For  a  distribution  that  can  be  described  by  Eq.  (1.  5),  the  total  number  of 
hydrometeors  between  a  minimum  diameter  (d)  and  a  maximum  (D^j)  is 

D„, 


NT=  i. 


Number  m 


N  r„ 
o  N 


Number  m 
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where  r^  is  the  truncation  ratio  of  the  number  of  hydrometeors  contained  within 
the  d  and  B^j  limits  to  the  total  N  as 


D 


M 


N 


I  N  dD 
Jd 
» 

J  N  dD 


or 


rN=e 


-dA  -DM  A 


(1.9) 


(1.10) 


The  hydrometeor  liquid  water  content  is  distributed  with  diameter  as 

MD  =  fx  10'3Pu)  NoD3e'AD  g  m”3  mm-1  .  (1.11) 


and  the  total  M  from  d  to  is 
DM 

M  =  J  Md  dD  g  m“3 
d 


or 


jX  10  3  p  N  rA.  „ 

o  M _ -3 

- - 3 -  g  m 


(1.12) 


(1.13) 


where  r^  is  the  truncation  ratio  of  the  liquid  water  content  contained  within  the 

limits  d  and  D„,  to  the  total  M  as 
M 


Jd  MDdD 

(1*  14) 

J  MD  dD 
o 


or 
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(1.  15) 
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i  |  e'dA  |(d  A  )3  +  3(d  A  )2  +  6dA+6] 


-D,.A 


e  M‘  (DmA)3  +  3  (DMM2+5DM«4  6l), 


The  radar- reflectivity  distributed  values  for  the  same  population  are 


„  n6  -nD  6  - 3  - 1 

Z„  =  N  D  e  mm  m  mm  , 
u  o 


(1.16) 


and  the  total  Z  within  the  limits  d  to  is 


D 
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z  =  r  ZD  dD 


6  -3 

mm  m 


(1. 17) 


or 


720  N  r2  6  _3 

Z  - - - -  mm  m 


(1.-18) 


where  r2  is  the  truncation  ratio  of  the  radar  reflectivity  contained  within  the  limits 


d  to  D  to  the  total  Z  as 

IVi 


D 


II 


/.  Zn<"> 


Vy  ~ 

00 


(1*  19) 


f  Z  dD 


rZ  ”  720 


|  e*dA  ((dA)6  +  6(dA)5  +  30(dA)4  +  120  (dA)3  +  360(dA)2  +  720dA+  720] 
-  e  °M  ((P^A)6  +  6(DmA)5  +  30(DMA)4  +  12C(DMA)3  +  360(DMA)2 


'll 


+  720(DMA)  +  720]  }  . 


M 
U.  20) 


When  Eqs.  (1. 13)  and  (1. 18)  are  solved  for  Nq  they  can  be  equated  as 


M  A*  Z  A' 

y  .  n-3  ""  720  r« 

w  X  10  Pu>  rn  2 


(1.21) 
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and  solved  for  A  as 


A  = 


M  rz  720 

- - - - 

Z  ir  X  10  p  r„. 

'w  M 


,1/3 


mm 


For  Pu}=  1  g  cm  #  this  reduces  to 


A  =  61.  2 


M  r. 


1/3 


mm 


-1 


M 


(1.22) 


(1.23) 


Since  both  r^  and  r^  incorporate  A  as  a  term  (Eqs.  (1. 15)  and  (1.  20)],  then  a  trial 
and  error  method  can  be  used  to  solve  Eq.  (1.  23)  whereby  the  value  of  A  is  adjusted 
until  both  sides  of  the  equation  are  equal.  Once  A  has  been  found,  Nq  can  be  deter¬ 
mined  using  either  Eq.  (1. 13)  or  (1. 18)  and  the  total  number  of  particles  may  be 
calculated  from  Eq.  (1. 8). 

If  the  M  and  Z  derived  from  the  truncated  precipitation  probe  measurements 
of  the  1-D  instrument  are  used  in  these  equations  with  the  appropriate  diameter 
limits,  where  DM  now  becomes  the  upper  truncation  diameter  DT»  then  the  dis¬ 
tribution  properties  of  the  -spectra  can  be  described  by  a  function  of  exponential 
type.  The  plotted  exponentials  of  Figures  1. 1  and  1. 2  (solid  lines)  were  calculated 
in  this  manner  where  d  =  0. 133  mm  and  DT  =  1. 496  mm.  It  is  apparent  from  these 
figures  that  the  solid  line  representing  the  exponential  number  distribution  estab¬ 
lishes  the  trend  of  the  spectrum  and  an  estimation  of  that  portion  of  the  population 
missed  because  of  instrument  truncation  can  be  made  by  extrapolating  the  line  to 
some  larger  diameter. 

Now  the  second  part  of  the  problem  becomes  paramount  and  that  is  the  estab¬ 
lishment  of  a  new  D„,  value  which  reflects  the  maximum  size  of  the  non-truncated 
spectra.  Fast  studies  of  observed  versus  exponential  distributions  conducted  at 
AFGL  have  shown  that  can  be  related  to  A  through  Eq.  (1. 6).  The  dimensionless 
constant  C  in  Eq.  (1.  6)  is  dependent  upon  crystal  type  and  density  and  remains 
relatively  constant  for  any  specific  snow  type  occurring  in  a  particular  storm.  The 
normal  range  of  C  is  between  9  and  12. 

When  the  maximum  particle  sizes  of  known  crystal  habit  that  are  present  in  any 
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particular  pass  are  measured  with  the  Aluminum  Foil  Hydrometeor  Sampler  and/or 
the  PMS  2-D,  they  can  be  converted  to  values  through  Eq.  (1.4)  which,  in  turn, 
can  be  used  in  Eq.  (1.6)  with  the  A  of  that  pass  to  define  C.  Thus,  the  relationship 

14.  Church,  J.F.,  Poes,  K.K.,  and  Spatola,  A.A.  (1975)  The  Continuous  Alumir 
num  Foil  Sampler;  Design  Operation,  Data  Analysis  "Procedures,  and 
Operating  Instructions,  Instrumentation  Papers  No.  235,  AFCRL-? R-75- 
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of  D  and  A  for  the  complete  pass  can  then  be  used  to  define  the  D  of  the  trun- 
m  m 

cated  spectra  for  any  portion  of  the  pass. 

The  maximum  physical  sizes  recorded  by  the  2-D  for  the  flight  passes  plotted 
in  Figures  1-  1  and  1.  2  were  approximately  8  and  6. 5  mm.  When  these  crystal 
sizes  were  substituted  in  Eq.  (1. 4)  and  the  parameters  of  £  and  0  were  given  the 

7 

values  of  0.4  and  0.782  which  are  designated  for  large  snow,  the  maximum  equivs 
lent  melted  diameters  of  2.03  and  1.  73  mm  resulted.  These  maximum  diameters 
were  then  used  in  Eq.  (1.  6)  with  the  A  values  from  the  passes  of  Figures  1. 1  and 
1.  2  to  determine  the  D  A  constants.  The  C's  for  both  passes  were  approximately 
9.  0. 

Since  N  and  A  are  known  and  a  IX,  has  been  established,  the  estimated  values 
o  M 

of  M  and  £  can  now  be  determined  by  finding  the  M  and  Z  contribution  from  the 
extrapolated  part  of  the  distribution  (dotted  line  in  Figures  1. 1  and  1.2)  and  adding 
them  to  the  measured  values  as 


ME  =  Mp  +  AM 


(1.24) 


ZE  =  Zp  +  AZ 


6  -3 

mm  m 


(1.25) 


where  M£  and  Z„  are  the  estimated  liquid  water  content  and  radar  reflectivity, 

Mp  and  Z^  are  the  values  derived  from  the  precipitation  probe,  measurements  and 

AM  and  AZ  are  the  contributions  fromthc  extrapolated  D_,  to  D,,  portion. 

I,  IV1 

The  equations  for  M„  and  Z„  thus  become 

it.  it. 


ME  ■  l  J,D  dD  +  X  JID  <*D  s 

d  u,.. 


(1.26) 


ZE=  l  ZDdD  +  /  ZDdD 
d  D.t> 


•  6  “3 

mm  m 


(1.27) 


Since  both  first  terms  are  known  quantities,  the  solving  of  these  equations  be¬ 
comes  a  two-part  process  where  Mp  and  Zp  are  used  in  Eq.  (1.  23Fto  determins 
the  A  of  the  exponential  function  which,  in  turn,  is  used  in  either  Eq.  (1.  13)  or 
(1. 18)  to  find  Nq.  These  parameters  are  then  incorporated  into  the  second  term  of 
Eq.  's  ( 1.  26)  and  ( 1.  27)  to  find  AM  and  AZ. 


Table  1.2  lists  the  Alp  and  Zp  that  were  derived  from  the  measurements  from 
the  precipitation  probe  for  the  two  sample  passes  compared  with  the  estimated 
Mg  and  Zp  from  the  extrapolated  spectra. 


Table  1. 2.  Comparison  of  Precipitation  Probe  Data  Mp,  Zp  With 
the  Estimated  Values  (Mp,  Zp)  Determined  from  the  Extrapolated 
Distributions 


MP 

-3 

g  m 

mf 

-3 

g  m 

FACTOR 

A1e/Mp 

Zp 

6  -3 

mm  m 

6E 

mm  m  -3 

FACTOR 

ZE/ZP 

Pass  3 

0.3824 

0.4144 

1.084 

544.  1 

845.5 

1.554 

Pass  5 

0.  3062 

0.3137 

1.  024 

336.4 

393.6 

1. 170 

1.3  CONCLUSIONS 

Although  the  amount  of  M  and  Z  missed  because  of  instrument  truncation  could 
vary  considerably  in  different  situations,  it  is  clear  from  these  results  that  it  is 
the  Z  parameter  which  shows  the  largest  degradation.  In  the  case  of  a  slightly 
truncated  distribution,  the  measured  M  may  be  acceptable  but  if  the  sampled 
hydrometeors  are  to  be  defined  in  a  relationship  of  aircraft  measured  Ai  and  2, 
then  new  estimated  values  of  Z  have  to  be  incorporated. 

The  reader  is  referred  to  Berthel3  for  a  more  detailed  discussion  on  how 
instrumentallv  truncated  aircraft  data  may  be  utilized  in  regression  analysis. 
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2.  Nominal  Equations  to  Determine  Liquid 
Water  Content  ( Mass)  From  Radar 
Reflectivity  fit  Kwqjalein 

Roberto.  Bertnet 


2.1  INTRODUCTION 

Knowledge  of  the  mass  of  liquid  water  CM)  encountered  Uy  a  missile  is. crucial 
in  the  evaluation  of  weather -related  missile  erosion.  In  tests  conducted  at  the 
Kvvajalcin  Missile  Range  (KMR),  the  M  profile  is  established  from  the  reflectivity 
measurements  (Z)  of  the  weather  radar  taken  along  the  trajectory  path  of  the  missile 
which  are  then  associated  with  the  aircraft-radar  relationships  of  M  vs  Z  as  deter¬ 
mined  from  that  specific  time  period.  These  relationships  are  developed  from 
simultaneous  radar  and  aircraft  measurements  of  cloud  and  precipitation  particles' 
at  various  altitude  levels  made  before  and  after  missile  flight.  The  number  densir  \ 
ties  measured  with  the  aircraft  instrumentation,  modified  by  other  meteorological 
information,  are  converted  into  "k"  factors  (defined  by  the  relation  k  =  M/^Z)  and 
are  then  correlated  with  the  radar  returns  to  give  power  function  relationships 
which  are  readily  converted  into  mass. 

Until  specific  aircraft -radar  equations  are  established,  any  determination  of  M 
from  radar  Z  is  necessarily  based  on  standard  relationships.  Up  to  the  time  of 
this  study,  nominal  equations  had  not  been  derived  from  the  KMR  area  and  the 
equations  being  used  as  stand  -ds  were  those  that  were  used  in  the  weather  erosion 
tests  conducted  at  Wallops  Island,  Virginia. 

15.  Plank,  V.G. ,  Berthel,  R.O. ,  and  Barnes,  A.  A.  (1980)  An  improved  method 
for  obtaining  the  water  content  values  of  ice  hydrometeors  from  aircraft 
and  radar  data,  J.  Appl.  Meteor.  19:1293-1299. 
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The  purpose  of  this  investigation  is  to  establish  standard  equations  peculiar  to 
the  climatological  conditions  of  the  Kwajalein  Atoll. 

2.2  DATA 

All  aircraft-radar  correlated  passes  that  were  taken  at  KMR  from  January  1977 
through  January  1979  were  considered  in  this  study.  The  original  data  tapes  from 
the  Learjet  (Aeromot  Inc..  Tulsa,  Oklahoma),  the  MC-I30E  (AFGL)  and  the  ALCOR 
and  TRADEX  radars  (from  MIT  Lincoln  Laboratory)  were  reprocessed  through  the 
AFGL  computers.  Each  separate  pass  was  then  reanalyzed  in  accordance  with  the 
following  criteria: 

Aircraft 

(a)  Instruments  operating  properly 

!b)  Particie  typing  (verification  of  original  type  or  change) 

(c)  Level  altitude  (aircraft  not  ascending  ot  descending; 

(d)  Pass  time  limits  (to  conform  to  periods  of  level  flight) 

Radar 

(a)  Instrument  operating  properly 

<b)  Correct  altitude,  range  and  elevation 

(c)  Reflectivity  did  not  include  returns  from  aircraft 

(d)  Sufficient  variability  in  Z  to  insure  a  correct  regression  analysis. 

The  few  passes  that  did  not  meet  these  criteria  were  discarded. 

Forty -three  passes  were  deemed  to  be  of  sufficient  quality  to  use  in  the  analysis: 

14  June  1977  6  passes  7  January  1978  4  passes 

6  August  1977  1  pass  23  June  1978  3  passes 

16  September  1977  4  passes  4  July  1978  17  passes 

4  November  1077  6  passes  26  June  1979  2  passes 

There  were  11  cases  of  Bullet-Rosettes,  21  cases  of  Small  Snow,  8  cases  of 
Large  Snow  and  3  cases  of  Rain.  Out  of  the  43  cases  used  in  this  analysis,  17  or 
39.5%  were  from  the  "heavy"  weather  test  of  4  July  1978,  the  remainder  being 
composed  of  "light"  weather  or  cloud  situations. 

2.3  DATA  ANALYSIS 

In  the  normal  analytical  procedure  used  at  AFGL,  the  multi-sample  spectral 
data,  acquired  by  the  Lear  and  MC-130E  along  particular  flight  paths  through  cloud 
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and  storm  situations,  are  used  to  obtain  equations  that  mathematically  define 
hydrometeor  environments.  This  is  done  by  the  conventional  method  of  cross  - 
plotting  the  logarithmic  values  of  M  and  Z  from  the  multiple,  individual  samples. 
The  least -square  analyses  of  this  field  of  data  points  determines  the  line  of  best 
fit  and  provides  the  coefficient  and  exponent  of  the  power  function  regression  equa¬ 
tion.  The  ultimate  goal  is  to  relate  ?.l  and  Z  in  a  regression  relationship  as 

M  »  a  Zb  (2.  1) 


so  that  the  liquid-water-content  of  a  specific  hydrometeor  region  can  be  calculated 
from  a  given  radar  reflectivity. 

Two  methods  have  been  employed  in  the  past  to  acquire  M-Z  information.  One 
is  the  regression  of  the  dependent  aircraft  derived  variables  of  versus  Z^,  and 

the  other  the  regression  of  the  indeoendent  variables  M  ,  versus  the  radar  measured 

„  3 

R’ 

Obtaining  the  M  and  Z  from  aircraft -measured  ice  hydrometeors  requires  the 
knowledge  of  the  sizes  of  water  drops  that  would  be  formed  if  the  ice  particles  were 
melted.  *  These  equivalent  melted  diameters  are  determined  by 

D  =  0  f  ^  (2.2) 


where  t  is  the  measured  physical  size  of  an  ice  particle  and  a  and  are  values 
assigned  to  a  particular  type  of  crystal. 

The  M  of  a  measured  distribution  is  found  by 

I  =  n 


)i, 


10 


N  a3  2 
1  i=  1 


32  .. 
i  **i 


g  m 


-3 


(2.3) 


where  P  is  the  density  of  water  (1  g  cm  )  and  y.  is  the  ratio  of  the  number  of 
particles  in  a  class  (i)  to  the  total  number  (Nj)  in  the  distribution. 

The  Z,  which  is  equivalent  to  the  reflectivity,  measured  by  a  radar,  from  a 
distribution  is  determined  by 


i  =  n 

■7  ~  m  -t®  T  t  Sd 

A  '  NTa  jT'j  £s  yi 


6  -3 

mm  m 


(2.4) 


These  equations  show  that  the  aircraft  spectral  M,  andZ  are  highly  dependent 
upon  the  a  and  .8  used  lo  convert  the  measured  physical  sizes  of  the  ice  particles 
into  equivalent  melted  diameters.  When  these  values  are  used  in  the  regression 
analysis  Eq.  (2. 11,  any  uncertainty  in  the  knowledge  of  crystal  type  is  reflected  in 
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both  Z  and  M.  Using  in  place  of  the  Z^  in  Eq.  (2. 1)  removes  these  particular 
uncertainties  from  one  parameter  of  the  analysis  since  Z^  does  not  require  a  prior 
knowledge  of  crystal  type  but  only  if  the  particles  are  water  or  ice. 

Past  studies  at  AFGL11'  have  shown  that  a  spectral  parameter  k,  which  is 


defined  as 


.  A  -3  -1.5 

k  -  . -  g  mm  m 

VZA 


when  used  in  the  regression  analysis  as 


will  minimize  the  uncertainty  effect  of  particle  typing.  The  related  M  of  the  sample 
is  then  found  bv 


b ,  0. 5 

c  a  7  1 

ai  R 


•12. 7) 


Each  correlated  pass  was  computer  processed  using  3 -sec  running  mean  aver¬ 
ages  to  produce  individual  k relationships  [Eq.  (2. 8)J.  The  cases  in  each  par- 
tide  type  category  were  then  combined  with  each  case  weighted  equally  to  form  a 
composite  equation  to  relate  51  and  {Eq.  <2.?|J. 

The  number  of  points  used  in  each  category  were: 


Bullet  -Rosettes 
Small  Snow 


Large  Snow 


2.4  RESULTS 


The  composite  equation  for  Bullet -Rosettes  is 


?.l  =  0.  04483  Z 


0.  4602 


with  a  standard  deviation  of  3. 13  dBM  (a  factor  of  2. 06.  ), 


Small  Scow: 


with  a  standard  deviation  of  2.  C"  dBM  (a  factor  of  2. 03. ), 
for  Large  Snow: 

M  =  0.02024  ZR°*  4406  (2.101 

with  a  standard  deviation  of  1.  *i2  dBM  (a  factor  of  1. 55.  ), 
for  Rain: 

M  =  0.005476  ZR°* 5,649  (2.  11) 

with  a  standard  deviation  of  2.79  dBM  (a  factor  of  1.90. ). 

As  previously  mentioned,  there  are  only  three  cases  of  rain  considered  in  this 
analysis.  The  procurement  of  more  correlated  data  would  be  expected  to  change 
this  equation  although  no  significant  differences  are  anticipated. 

All  eight  cases  of  large  snow  were  obtained  on  the  same  date,  4  July  1978, 
and  seven  of  these  cases  exhibit  some  degree  of  instrument  truncation.  When  the 
extrapolation  procedure  outlined  in  Section  l  is  used  on  these  data,  the  composite 
equation  for  large  snow  [Eq.  (2. 10)]  becomes 

M  =  0.  02463  2o0*  3890  .  (2-  12) 

Comparison  of  M  calculated  from  the  equation  derived  from  truncated  data  [Eq. 

(2.  10)]  and  from  the  extrapolated  Spectra  {Eq.  (2. 12)]  show  that  the  noncorrected 
equation  gives  higher  values.  For  a  Z  =  300,  the  increase  is  9.7%,  and  when 
Z  =  700,  the  increase  is  14.  6%. 

It  must  be  emphasized  that  these  nominal  equations  are  averages  that  were 
obtained  using  all  available  data.  The  standard  deviations  (and  factors)  include  67% 
of  these  data  thus,  33%  vf  the  data  He  outside  of  the  Umits  and  therefore,  in  some 
cases,  the  data  points  deviate  considerably  from  the  average. 

The  computer  processing  of  radar  reflectivity  using  these  nominal  equations 
will  give  quick  and  convenient  estimations  of  liquid-water-contents  for  use  in  the 
predictions  of  specific  weather  criteria  when  the  sampling  aircraft  is  unable  to  be 
on  station  or  when  only  crystal  habit  information  is  available  from  aircraft.  Use 
of  these  equations  are  not  intended  to  be  replacements  for  those  individualM-Z 
relationships  that  are  derived  from  simultaneous  aircraft  and  radar  measurement 
for  specific  time  periods. 
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3.  Variation  of  Mk  With  Altitude 
and  With  Temperature 

Y.tZMg 


3.1  INTRODUCTION 

The  problem  of  combining  aircraft  and  radar  data  Jed  to  the  introduction  of  the 
factor  k.  where  k  *  Miy  2..  In  the  coarse  of  data  collection  at  Watimss  Island  awl 
at  Kwajalein,  It  was  noted  that  k  varies  with  particle  type;  average  values  Increase 
as  we  progress  from  rain,  to  large  scow,  to  small  snow,  and  finally  to  ballet 
rosettes.  Because  particle  type  is  to  a  large  extent  a  function  of  temperature  aM 
therefore  altitude,  it  was  decided  to  search  for  some  relation  between  k  and  alti¬ 
tude,  or  between  k  and  temperature.  "  * 

This  section  presents  the  results  of  analyzing  data  from  Wallop  Island  and 
Kwaialein.  and  expressing  k  as  a  function,  first  of  altitude,  then  of  tern  pc  r  at  ore. 

In  order  to  simplify  the  presentation,  the  parameter  Ilk  ®  1000  k  Is  used.  Mk  ranges 
from  approximately  1  to  ISO. 

3.2  BATA 


The  Mk  values  prnsented  here  are  pass  averages  of  3 -sea  running  means  ob¬ 
tained  during  aircraft  and  radar  correlation  runs.  The  Wallops  correlation  runs 
were  made  on  eight  different  days  during  the  period  from  January  through  April 
1911,  and  the  Kwajalein  runs  were  made  on  seven  different  days  during  the  iS-month 
period  starting  June  1871. 
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A  number  of  correlation  runs  were  deleted  from  this  study  to  avoid  misinter¬ 
pretation.  Those  deleted  were  as  follows: 

11)  Runs  with  missing  or  extremely  small  values  of  31.  The  total  aircraft 
M^,  which  is  the  stun  of  the  cloud  31^.  and  the  precipitation  Sip,  and  winch  is 
directly  related  to  31k,  is  not  reliable  in  such  eases,  particularly  In  the  regime 
where  the  contribution  of  cloud  to  the  total  M  is  expected  to  be  large. 

(2)  Runs  in  which  the  comparison  of  the  arithmetic  mean  with  the  geometric 
mean  of  the  aircraft  data  showed  large  discrepancies.  In  such  cases  the  variability 
of  the  data  was  too  great  to  obtain  meaningful  averages. 

A  total  of  1 10  Mk  values  were  thus  available  for  this  study,  with  65  from 
Wallops  and  45  from  Kwajalein. 

3  3  RESULTS 

The  variation  of  Mk  with  aircraft  altitude  for  Wallops  is  shown  fit  Figure  3. 1  and 
for  Kwajatein  in  Figure  3.2.  The  variation  of  Mk  with  aircraft  measured  tempera¬ 
ture  for  Wallops  is  shown  in  Figure  3.3  and  for  Kwajalein  in  Figure  3.4.  The 
particle  types  are  indicated  In  all  four  figures.  The  figures  show  a  much  greater 
range  diuiuue  and  temperature  at  KwajateUV  than  at  W3ii"p«  Island;  The  data 
points  scatter  widely  about  the  regression  lines,  but  the  trend  is  far  Mk  to  increase 
with  increasing  altitude  and  to  decrease  with  increasing  temperature.  The  re¬ 
gression  equations  obtained  are  listed  in  Table  3.  1.  In  the  equations  B  is  the  height 
In  kilometers  and  T  Is  the  temperature  in  degrees  Celsius.  Along  with  the  equa¬ 
tions  are  tabulated  the  correlation  coefficients  and  the  rms  values. 

Because  of  the  different  temperature  and  altitude  relationships  sound  In  the  two 
markedly  different  climatic  regions  of  Wallops  Island  <75aW.  3SSN)  ami  Kwajalein 
i!iS33E.  9*N).  signifleant  differences  were  expected,  particularly  between  the  height 
equations.  However,  the  results  showed  surprising  similarity. 

Data  were  obtained  up  to  greater  heights  and  there  was  a  larger  range  of  heights 
sampled  at  Kwajalein.  Nevertheless,  the  differences  in  the  Mk  -  Height  relation  r 
ships  at  Wallops  and  Kwajalein  are  meteorologically  significant.  As  can  be  seen  by 
comparing  Figures  3. 1  and  3. 2,  these  differences  were  larger  at  the  higher  heights 
and  undoubtedly  were  related  to  the  more  extensive  depths  of  the  storms  at  Kwajalein. 
Wallow  Island  storms  topped  out  between  1  and  10  km  since  the  storms  sampled 
were  limited  to  winter  and  early  spring,  when  the  tropopause  was  generally  below 
10  kilometers.  With  Kwajalein  so  close  to  the  equator  there  was  less  variation  of 
tropopause  height,  18  to  17  kilometers. 

The  spread  of  Mk  values  versts  temperature  for  the  various  crystal  habits  was 
less  at  Wallops  Island  than  it  was  at  Kwajalein,  In  addition,  H  should  also  be  noted 
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that  the  correlations  were  smaller  and  the  rms  values  were  larger  at  Kwajalein  as 
seen  in  Table  3. 1,  This  can  be  attributed  to  the  greater  extent  and  more  convective 
nature  of  the  storms  at  Kwajalein. 

Figures  3.  3  and  3. 4  clearly  show  different  relations  between  Mk  and  tempera¬ 
ture  at  the  two  locations  thus  indicating  that  different  relationships  should  be  used 
for  tropical  and  extratropical  storm  systems. 

The  above  results  show  that  reasonable  estimates  of  Mk  as  a  function  of  height 
or  temperature  can  be  made  for  either  Wallops  Island  or  Kwajalein.  Whether  the 
results  are  accurate  enough  depends  on  the  accuracy  requirements  of  the  particular 
experiment. 


Table  3.  1.  Mk  as  a  Function  of  Height  or  Temperature  at 
Wallops  Island  and  Kwajalein 


Location 

Regression  Equations 

Correlations 

rms  . 

Wallops 

Mk  =  4.  699  e°* 284  H 

0.833 

0.433 

Kwajalein 

Mk  =4.919  e0,205  H 

0.772 

0.487 

Wallops 

Mk  =  12.217  e~°'  0474  T 

-0.868 

0.388 

Kwajalein 


-0.  0291  T 


Mk  *  14.  264  e 


-0.739 


0.516 


TEMPERATURE  (*C) 


4.  Predicting  Trajectory  M  Values  at  Wallops  Island 


R.  M.  Dyer 

4.1  INTRODUCTION 

Because  it  is  not  possible  to  make  aircraft  measurements  along  the  trajectory 
at  the  time  of  missile  reentry,  a  relation  between  radar  Z  and  trajectory  M  has  to 
be  assumed.  At  first,  this  was  a  power  relation  of  the  form  M  =  aZ*\  with  a  and  b 
obtained  from  results  published  by  observers  at  locations  other  than  either  Wallops 
Island  or  Kwajalein.  As  data  were  acquired  at  Wallops  Island,  especially  tailored 
M-Z  relations  were  derived  for  each  crystal  type  shortly  before  and  shortly  after 
each  reentry.  Assuming  no  temporal  change,  these  relations  were  then  applied  to 
the  reentry  time. 

From  the  initial  Wallops  Z_  data  using  the  literature  M-Z's,  it  was  possible 

K  16 

to  derive  a  climatology  of  Wallops  Island  storms.  Climatological  values  of 

k  (k  )  can  be  obtained  using  these  data  and,  once  k  has  been  derived  for  each 
c  c  0  5 

particle  regime,  a  climatological  -  Z^  relation  of  the  form  s  Z^  '  can 

then  be  applied  to  future  situations  when  radar  data  are  available. 

Several  questions  then  arose.  Can  measured  aircraft  (M^)  and  radar  (Zp)  rela 

tiottr.  be  applied  to  a  meteorological  situation  several  hours  before  or  after  the 

measurements  are  made?  How  does  this  differ  from  a  M-Z  relation  obtained  from 
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previous  measurements  reported  in  the  literature?  Does  the  accuracy  improve  if 
climatological  values  of  It  are  used? 

Finally,  are  the  ranges  of  M  encountered  at  a  given  altitude  or  temperature, 
or  associated  with  a  given  particle  type,  sufficiently  narrow,  and  is  there  an 
adequate  amount  of  data,  that  an  appropriate  climatological  value  of  Mc  could  be 
substituted  for  radar  or  aircraft  measurements?  The  latter  question  applies  to 
missile  reentry  over  areas  for  which  it  is  not  possible  to  have  either  radar  or 
aircraft  measurements  on  the  spot. 

The  analysis  reported  in  the  following  pages  is  an  attempt  to  answer  these 
questions  for  Wallops  Island.  The  methodology  can  be  applied  to  other  locations 
as  well. 

4.2  THE  DATA 

During  the  first  four  months  of  1977,  special  efforts  were  made  at  Wallops 
Island,  Virginia,  to  make  correlated  aircraft  runs;  that  is,  to  measure  the  cloud 
physics  parameters  along  a  constant  altitude  path  at  the  same  time  that  a  ground- 
based  radar  was  measuring  the  radar  reflectivity  along  the  path.  Each  run  pro¬ 
duced  an  aircraft -radar  relationship  of  MA  -  2^. 

Several  successive  aircraft -radar  correlation  runs  were  made  at  virtually  the 
same  altitude,  at  intervals  ranging  between  a  few  minutes  and  three  hours.  Each 
data  point  consisted  of  a  4 -sec  sample  of  the  hydrometeor  spectrum  obtained  from 
the  PMS  1-D  probes  aboard  the  aircraft,  and  the  measured  ground-based  radar 
return  from  the  volume  of  space  one  range  gate  short  of  the  gate  containing  the 
aircraft  echo.  A  correlation  run  was  useable  for  the  study  reported  here  if  and 
only  if  more  than  one  pass  were  made  on  the  same  day  within  the  same  altitude 
interval,  and  with  the  same  predominant  particle  type. 

Useable  correlation  runs  were  made  on  five  separate  days  between  9  January 

and  4  April,  1977.  Altogether,  there  were  38  constant  altitude  passes  for  which 

both  aircraft  and  radar  data  were  available.  The  passes,  divided  into  four  broad 

categories  of  particle  type,  are  tabulated  in  Table  4.  1.  For  the  purpose  of  this 

analysis,  they  were  grouped  into  pairs,  A  pair  consists  of  two  correlation  passes 

made  on  the  same  day,  at  the  same  altitude  and  with  the  same  observed  particle 

type.  The  -  7, ^  relationship  derived  from  each  run  was  then  applied  to  the 

pass -average  ZR  of  its  pair  partner,  which  gave  a  derived  M  and  then  was  compared 

with  the  partner's  observed  AT.  Using  only  the  criteria  of  same  day,  same  altitude, 

and  same  particle  type,  it  was  possible  to  obtain  42  pairs  of  correlated  pass  data 

that  could  be  used  to  test  the  predictive  power  of  the  M  „  -  Z~  relations  derived 

a  tt 

from  the  combination  of  aircraft  and  radar  data.  However,  only  18  of  these  were 


completely  independent,  in  the  sense  that  each  pass  was  used  only  once.  As  a 

first  exercise,  only  independent  pairs  were  used.  Each  pass  was  paired  with  the 

pass  closest  to  it  in  time,  adhering  to  the  same  altitude,  same  particle  tvpe 
requirement. 


Table  4.  1.  Summary  of  Correlated  Aircraft  Radar  Data 


Particle  Type 

Altitude 

Passes 

Pairs 

Independent  Pairs 

Rain 

0.  5  km 

4 

6 

2 

2.  1  km 

4 

6 

2 

3.  1  km 

2 

1 

) 

TO 

TO 

5 

Large  Snow 

3.  8  km 

3 

3 

1 

3.  9  km 

2 

1 

1 

4.8  km 

4 

6 

2 

S' 

TO 

¥ 

Small  Snow 

5.  8  km 

4 

6 

2 

6.  0  km 

2 

1 

1 

6.  3  km 

3 

3 

1 

7.  0  km 

4 

6 

2 

TO 

TO 

F 

Bullet  Rosettes 

7.  9  km 

2 

1 

1 

8.  5  km 

2 

1- 

1 

9.  3  km 

2 

1 

1 

F 

3 

3 

Totals 

38 

42 

18 

4.3  ANALYSIS  PROCEDURE 

Five  methods  of  estimating  liquid  Wi..er  content  at  a  point  were  evaluated  in 
this  study. 

(1)  In  the  absence  of  any  measurements,  whether  radar  or  aircraft,  it  is 
necessary  to  rely  on  climatological  data.  Because  ZR  measurements  were  made 
between  1971  and  1975  at  Wallops  Island,  there  is  an  adequate  data  base  for  deriving 
an  average  TO,  as  a  function  of  hydrometeor  type  and  synoptic  situation,  using 
literature  M-Z  relations.  These  averages,  published  by  Berthel, iS  were  compared 
with  the  pass  average  M  obtained  from  the  1977  correlated  aircraft -radar  data. 
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Since  a\\  »*•"  ccrreialiun  runs  were  made  in  stratiform  storm  situations,  the 
climatological  values  were  taken  from  Table  A1  for  stratiform  storms  from 
Berthel's  report.  They  are  listed  in  Table  4.  2. 


Table  4.  2  Climatological  Average  Liquid  Water  Content  Values  at 
Wallops  Island  Based  on  Zp  and  Literature  M-Z  Relations  for 
Stratiform  Situations.  Taken  from  Table  A1  of  Berthel  (1976) 


Hydrometeor  Type 

_3 

Mean  (gm  m  ) 

Standard  Deviation 
(gm  m"3)  SD/Mean  (%) 

Bullet  Rosettes 

0.  087 

0.073 

83.9 

Small  Snow 

0.  116 

0.  052 

44.8 

Large  Snow 

0.  257 

0.  124 

48.2 

Rain 

0.  161 

0.  095 

59.0 

The  large  standard  deviations  indicate  that  the  climatological  means  for  strati¬ 
form  systems  will  not  be  accurate  indications  of  liquid-water-content,  and  cer¬ 
tainly  are  not  desirable  if  either  radar  or  aircraft  data  are  available.  However,  if 
Wallops  Island  were  in  a  location  without  satellite  or  conventional  weather  data, 
climatological  data  would  be  the  only  available  estimator.  It  is  useful,  therefore, 
to  determine  what  errors  can  arise  under  such  circumstances. 

(2)  In  the  absence  of  climatological  data,  and  with  only  a  ground-based  radar 
available,  recourse  must  be  made  to  M  -Z  relations  of  the  type  presented  by 
Plank.  ^  These  were  obtained  from  the  published  literature,  without  reference  to 
synoptic  situation  or  geographic  location.  The  derivation  of  the  literature  M-Z 
relations  is  shown  at  the  extreme  right  of  Figure  4.  1.  The  absence  in  tl.e  diagram 
of  intermediate  calculations  and  derived  quantities  between  the  selection  of  crystal 
type  and  the  ultimate  statistic  (the  Mp  -  ZR  relation)  does  not  indicate  that  these 
computations,  with  their  attendant  assumptions  and  sources  of  error,  were  not  re¬ 
quired,  but  only  that  we  have  not  tabulated  the  different  assumptions  of  the  various 
authors  who  derived  the  equations  originally.  The  relations  used  arc  listed  in 
Table  4.3. 

(3)  If  a  ground-based  radar  is  available,  and  a  sufficient  body  of  data  has  been 
collected  at  the  same  site  in  the  past,  it  is  possible  to  derive  an  M-Z  relation  for 
the  storm  type  and  geographic  location  of  interest.  The  steps  used  in  this  deriva¬ 
tion  are  shown  in  Figure  4,  1,  under  the  heading  "Previous  Measurements  (Clima¬ 
tology)".  Using  an  appropriate  M-Z  relation  from  the  literature,  each  radar  return 

Z  is  converted  to  a  measure  of  M.  The  climatological  k  should  be  determined 

_  0  5  ^ 

from  the  individual  ks,  that  is,  kc  3  M/Z  .  The  original  data  from  Wallops  Island 


were  no  longer  available  so  the  climatological  k  values  were  obtained  from  the 


average  M  (51  )  and  average  Z  (Z  }  values10  using  1?  =  M_/Z~ 


Using  these 


climatological  k  values  an  estimate  of  M  can  then  be  obtained  from  the  observed 
c  0  5 

Z„  values  and  the  relation  M  =  It  Z-  *  . 

R  c  R 

The  climatological  kc  values  used  here  were  derived  from  the  Berthel  data. 
For  bullet  rosettes,  5Ht  (k  X  1000)  is  65;  in  small  snow  ME  is  31;  in  large  snow 

C  C  'C 

19;  and  in  rain,  4. 


Table  4*  3.  M  -Z  Relations  Published  in  the 
Literature  and  Applied  to  Wallops  Measurements 


Bullet  Rosettes 

M  =  0.38 

20. 529 

Small  Snow 

M  =  0.0145 

z0.  538 

Large  Snow 

M  =  0. 00495 

20. 596 

Rain 

M  =  0.00314 

20.  576 

FROM  MEASUREMENTS  FOR  EACH  PASS 


PREVIOUS  MEASUREMENTS 
(CLIMATOLOGY! 


(LITERATURE! 


(CMJSERVBDJ  P*rtkkSp*ctmm  |  Rad* R«tom  foetid* Sftfcttim  |  1  RtdmHtwn  I  I 


DERIVED  I  Akers*  Z  {  |  AirtxaltM 


INTERMEDIATE 


SECOND  GENERATION 
DERIVED 


ULTIMATE 

STATISTIC 


Avcrsft  Zp 


Oimtolo^kM  OimMoteiiCTl 

iof  Mdt  crymi  Zr  -  Mft 

type  tesacbtrp** 

iyp* 


lot  Mdi  crystal 


Figure  4. 1.  Flow  Diagram  Showing  the  Methods  of  Deriving  M 


(4)  The  ideal  situation  is  to  have  both  aircraft  and  radar  data  available  at  the 
time  and  place  of  interest.  In  practice,  it  is  not  possible  to  make  aircraft 


measurements  along  the  reentry  trajectory  at  the  time  of  reentry.  The  next  best 
thing  is  to  make  measurements  as  close  to  the  reentry  time  as  possible,  then  to 
apply  the  M  -Z. relations  derived  during  the  calibration  pass  to  radar  measure¬ 
ments  taken  at  the  reentry  time.  This  was  done  for  each  of  the  passes  listed  in 
Table  4.  1.  The  power  relation  MA  =  aZ^  was  determined  for  each  pass,  and  this 
relation  was  then  applied  to  the  average  ZR  of  the  pass  which  was  paired  with  the 
calibration  pass. 

(5)  The  fifth  estimate  of  trajectory  M  is  the  pass  average  obtained  from 
the  regression  formula  derived  from  the  aircraft  and  radar  data  measured  during 
the  pass  itself  and  reapplied  to  Z^  values.  It  was  assumed  for  the  purpose  of  this 
study  that  these  data  represent  the  conditions  prevailing  during  the  missile  reentry. 
In  practice,  it  is  not  feasible  to  take  aircraft  measurements  during  reentry,  and  the 
closest  measurements  would  have  to  be  those  of  case  4.  However,  it  was  necessary 
to  establish  the  most  accurate,  or  reference,  value  of  M,  and  that  is  why  this  fifth 
method  is  included.  Even  this  reference  value  has  an  error  associated  with  it, 
namely,  that  resulting  from  applying  an  average  M-Z  relation  obtained  along  a 
constant  altitude  pass  to  any  point  on  that  pass.  This  is  the  minimum  error  of  any 
point  estimate,  and  results  from  the  variability  in  M  along  the  path.  It  is  expressed 
in  Table  4.5  as  the  standard  error  of  estimate  of  M  from  the  pass  derived  M  -Z 
relation. 


4.4  RESULTS 

Details  of  each  pass  and  the  results  of  the  computations  are  shown  'n  '•abies 

4.4a  through  4. 4d.  The  errors  (or  more  accurately,  the  uncertainties)  when  using 

each  method  of  prediction  are  expressed  as  percent  deviations  from  the  mean  M_, 

*  K 

The  standard  error  of  estimate  is  assumed  to  be  an  expression  of  the  minimum 
possible  uncertainty.  Whenever  a  predictive  method  produces  an  apparently  lower 
error  (those  values  in  parenthesis),  it  is  adjudged  to  be  as  accurate  as  meteorologi¬ 
cal  variability  will  allow.  Only  the  results  of  using  the  pass  average  next  closest 
in  time  are  listed  in  the  table. 
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Av«rn«  gundurd  Error  «  18.  BH 
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It  is  immediately  obvious  that  for  stratiform  storms  at  Wallops  Island  the 
competition  for  most  accurate  method  of  prediction  lies  between  using  correlated 
aircraft  and  radar  data  from  a  same-day  constant  altitude  pass  (ease  4),  and  using 
only  radar  data  and  a  climatological  value  determined  by  the  predominant  crystal 
type.  In  Table  4. 5,  the  average  improvement  in  accuracy  using  correlated  aircraft- 
radar  data  rather  than  climatological  !?,  values  is  expressed  as  the  average  change 
in  percent  error.  Negative  values  indicate  that  for  Wallops  Island  the  climatological 
method  vielded  better  correspondence  than  uie  pass  average  M-Z  or  Sc-Z 

relations.  Only  outlet  rosette  situations  <6  passes)  resulted  in  any  improvement 
when  aircraft  data  were  used.  These  results  were  not  changed  when  all  aircraft 
data  for  a  given  day  were  combined  to  make  the  prediction.  Nor  was  ther"  any 
improvement  when  the  appropriate  k  was  selected  from  F*jiii>es  3. 1  and  3. 3: 
according  to  the  altitude  or  ieniper—iut  e  at  which  the  pass  was  made.  However, 
using  the  average  k  measured  during  the  closest  correlated  pass  did  improve  the 
prediction  in  several  instances. 


Table  4.  S.  Changes  in  Percent  Error  of  M:  Climatological  k  Method  vs  Pass 

Regression  Method  c 


All  Passes 

Independent  pairs 
Only 

Correlation  With  Time 
From  Prediction 

Rain 

-4.5 

-6.S 

0.06 

Large  Snow 

-21.3 

-:  X  l 

-0..17 

Small  Snow 

-2.4 

-13.0 

-0.  50 

Bullet  Rosettes 

4-197. 

+  197. 

-0. 15 

Except  for  some  insignificant  changes  in  the  magnitude  of  the  errors,  the  re¬ 
sults  were  the  same,  whether  the  "average  error"  was  denned  as  the  average  of 
the  errors  in  the  individual  point  estimates  as  shown  in  Table  4, 5,  or  the  error  In 
the  estimate  of  the  average  M  along  the  pass.  In  either  case,  the  ranking  of  the 
estimates  remained  unchanged.  These  rankings  are  shown  In  Table  4.  6.  A  ranking 
of  i  was  given  to  Use  method  of  estimate  with  the  smallest  average  error.  Again, 
parenthesis  enclose. those  methods  of  estimate  which  yielded  an  average  error 
tower  titan  the  standard  error. 
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Table  4.  6a.  Relative  Rankings  of  Methods  of  Estimating  M  for  Rain 


Methods  of  Estimate 

Pass  Number 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Climatology 

5 

5 

4 

5 

5 

5 

5 

5 

(1) 

5 

Climatological  k^. 

(3) 

(2) 

(3) 

(3) 

(1) 

(i) 

(3) 

(1) 

(2) 

4 

Literature  M  =  aZ 

4 

4 

5 

4 

4 

3 

(4) 

4 

5 

3 

Nearest  Pass  M  =  aZb 

(2) 

(1) 

1 

(2) 

2 

4 

(2) 

3 

3 

2 

Nearest  Pass  k 

(1) 

(3) 

<1> 

(1) 

3 

(2) 

(1) 

(2) 

4 

1 

Table  4.  6b.  Relative  Rankings  of  Methods  of  Estimating  M  for  Large 
Snow 


Methods  of  Estimate 

Pass  Number 

1 

2 

3 

4 

5 

6 

7 

8 

9 

Climatology 

4 

3 

5 

2 

2 

1 

2 

3 

5 

Climatological  kc 

3 

4 

3 

1 

1 

3 

(1) 

2 

(2) 

Literature  M  =  aZb 

5 

5 

4 

3 

5 

4 

4 

5 

4 

Nearest  Pass  M  =  aZ 

(1) 

(1) 

2 

4 

4 

5 

5 

4 

(1) 

Nearest  Pass  k 

(2) 

(2) 

(i) 

5 

3 

2 

3 

1 

3 

Table  4.  6c.  Relative  Rankings  of  Methods  of  Estimating  M  for  Small  Snow 


Methods  of  Estimate 

t 

2 

3 

4 

5  6 

7 

8 

9 

10 

11 

12 

13 

Climatology 

0 

3 

2 

5 

5  4 

(4) 

4 

2 

(5) 

5 

5 

5 

Climatological  k^. 

(2) 

2 

3 

2 

(1)  3 

(3) 

2 

3 

(3) 

3 

4 

3 

Literature  M=  aZ*3 

4 

5 

4 

4 

4  5 

(5) 

5 

1 

(4) 

4 

(3) 

2 

Nearest  Pass  M=  aZ*3 

(3) 

4 

5 

3 

3  2 

(2) 

3 

4 

(2) 

(2) 

(1) 

4 

Nearest  Pass  k 

(1) 

(1) 

1 

1 

(2)  1 

(1) 

1 

5 

(1) 

(1) 

(2) 

(1) 

Table  4.  6d.  Relative  Rankings  of  Methods  of  Estimating  M  for 
Bullet  Rosettes 


Methods  of  Estimate 

Pass  Number 

1 

2 

3 

4 

5 

6 

Climatology 

4 

3 

5 

5 

3 

3 

Climatological  k^ 

5 

5 

4 

3 

5 

5 

Literature  M  =  aZb 

3 

4 

1 

4 

1 

4 

Nearest  Pass  M  =  aZb 

(2) 

(1) 

2 

(1) 

4 

(2) 

Nearest  Pass  K 

(1) 

(2) 

3 

(2) 

2 

(1) 

The  elapsed  time  between  the  pass  during  which  the  M  -  Z  relation  was  derived, 
and  the  pass  for  which  the  relation  was  applied  to  estimate  M,  varied  between 
three  minutes  and  three  hours,  two  minutes.  Before  this  study  was  undertaken, 
it  was  expected  that  the  rc curacy  of  the  estimate  would  be  highest  for  the  shortest 
elapsed  time  between  passes,  and  that  there  would  be  a  systematic  decrease  in 
accuracy  as  the  time  between  passes  increased.  This  was  not  the  ease.  As  the 
figures  in  the  last  column  of  Table  4,  5  show,  the  correlation  coefficient  between 
accuracy  and  time  between  the  derivation  of  the  M  -  Z  relation  and  the  prediction 
is  not  significant.  Indeed,  except  for  the  rain  cases,  this  coefficient  was  negative, 
indicating  that  if  the  correlation  were  significant  the  M  -  Z  relations  derived  from 
measurements  taken  only  a  few  minutes  before  the  prediction  would  be  less 
accurate  than  those  taken  hours  before!  This  conclusion  runs  counter  to  our  in¬ 
tuition.  Additional  analyses  using  independent  data  are  required  to  further  investi¬ 
gate  this  point  and,  at  present,  these  data  are  not  available.  It  should  be  restated 
here  that  all  of  these  correlation  runs  were  made  in  stratiform,  moderate  to 
heavy  intensity  storm  situations  during  winter  months  at  Wallops  Island,  Virginia, 
Attempts  to  characterize  predictability  based  on  synoptic  type  must  await  further 
data. 

Another  approach  would  be  to  use  M  -  Z  relations  derived  from  passes  which 
showed  the  same  characteristics  as  the  pass  which  is  to  be  predicted.  Because 
the  radar-measured  Z^  is  the  one  parameter  obtainable  in  real-time  during  reentry 
it  was  decided  to  test  the  hypothesis  that  the  M  -  Z  relation  obtained  during  an  air¬ 
craft  pass  could  be  applied  to  any  other  pass  with  the  same  hydrometeor  type  and 
approximately  the  same  average  Z^.  Data  to  test  this  hypothesis  proved  to  be 
very  sparse.  There  were  no  correlation  passes  in  bullet  rosettes  for  which  the 
average  Z  's  were  comparable.  Three  passes  in  small  snow  at  altitudes  between 
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4,7  and  6,4  km  had  average  Z„'s  between  25.  6  and  29,  2,  and  two  passes  in  large 

K  g  g 

snow  had  altitudes  and  ZB's  of  3.  9  km  and  149  mm  m  and  4.  8  km  and  155  mm 
-3  “ 

m  respectively.  The  power  relations  obtained  during  these  passes  were  used  to 
estimate  the  liquid  water  content  for  all  other  passes  having  approximately  the  same 
average  Z^.  In  the  case  of  rain,  there  were  two  passes  taken,  1  hr  10  min  apart, 
for  which  the  altitudes  differed  widely,  but  the  Z^'s  were  approximately  the  same. 

A  comparison  was  made  between  these  two  passes  to  determine  whether,  in  this 
case,  more  accurate  predictions  could  be  made  using  calibrations  obtained  at 
similar  Z_’s  rather  than  for  passes  close  to  the  time  and  altitude  for  which  the 

IT 

prediction  is  to  be  made. 

The  results  are  shown  in  Table  4.7.  For  the  rain  case  the  nearest  pass  was 
by  far  the  better  predictor.  In  the  case  of  large  snow,  using  the  same  Z^  pass 
data  with  the  closest  Z^  (obtained  more  than  a  month  previously)  was  preferable  to 
either  the  nearest  pass  or  the  climatological  It  calibration.  There  was  little  differ¬ 
ence  between  the  two  results  for  two  of  the  three  small  snow  passes.  For  the  third 
one,  a  highly  variable  case,  the  nearest  pass  calibrator  was  far  superior.  Though 
far  from  conclusive,  the  results  indicate  that  further  comparison  of  correlation 
runs  at  similar  altitudes  and  with  comparable  ZR's  should  be  undertaken  to  confirm 
these  conclusions  and  to  obtain  more  definitive  results. 

The  average  absolute  error  in  the  estimate  of  M  was  the  same  for  both  the 

_  -3 

climatological  k  and  the  nearest  pass  estimate,  0.  05  g  m  .  The  standard  devia- 

c  —  -3-3 

tion  was  higher  for  the  climatological  kc  predictions,  0.  10  g  m  vs  0.08  g  m 

for  the  nearest  pass  predictions.  In  general,  there  was  a  tendency  for  the  percent 

variability  along  a  pass  to  be  larger  for  the  smaller  average  M’s,  thus  pushing  the 

absolute  errors  toward  the  central  value  of  M.  However,  the  correlation  between 

percent  error  and  magnitude  of  M  (-0.  31  for  the  combined  data)  was  not  significant. 

4.5  CONCLUSIONS 

With  the  exception  of  the  b'Tlet  rosette  case,  there  now  seems  to  be  a  large 
enough  data  base  for  moderate  t**-  '  eavy  stratiform  storms  at  Wallops  Island  to 
permit  prediction  of  pass  average  liquid-water-content  using  the  radar  alone  and 
a  climatological  relation  br'ween  Z  and  M.  For  the  stratiform  storms  the  most 

stable  method  appears  to  be  ,o  use  a  climatological  value  F  ,  based  on  crystal 

0  5  c 

type,  and  the  relation  M  =  k  Z  . 


44 


It  cannot  be  emphasized  too  strongly  that  the  analysis  reported  here  does  not 

address  the  question  of  accuracy  in  measuring  M,  2  or  k.  That  subject  has  been 
17 

treated  elsewhere.  This  analysis  was  confined  to  a  comparison  of  different 
methods  of  estimating  M  along  a  trajectory,  assuming  that  all  input  parameters 
were  completely  error-fee. 

Under  ideal  conditions  (errorless  aircraft  data,  perfect  radar  measurements 
and  a  Z-M  relation  which  remains  constant  over  time  and  space),  there  would  still 
be  an  error  in  the  estimate  of  liquid  water  content  at  a  point  along  the  path,  due  to 
scatter  about  the  Z-M  regression  line.  The  deviations  of  a  and  b  result  in  errors 
in  the  estimate  of  M.  For  the  Wallops  island  data  used  here,  the  average  standard 
error  in  M  was  ±  IB.  8%  in  rain,  ±  10.  4?o  in  large  snow,  ±  23.  OTo  in  small  snow, 
and  ±  14.  1%  in  bullet  rosettes.  This  scatter  is  due  to  the  meteorological  variability 
inherent  in  even  stratiform  situations,  and  the  figures  cited  here  may  be  taken  as 
the  normal  atmosphere  variability  in  any  estimate  of  point  values  of  M. 


17.  Crane,  R.  K.  (1978)  Evaluation  of  Uncertainties  in  the  Estimation  of 
Hydrometeor  Mass  Concentrations  Using  Spandar  Data  and  Aircraft 
Measurements,  Sci.  Rep.  No.  1.  Contract  F19628-7(J-C-()069.  Environ  - 
mental  Research  and  Technology,  Inc.,  also  AFGL-TR -78-0118,  AD  A059223. 
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5.  Conclusions  and  Recommendations 


The  first  three  sections  of  this  report,  though  approaching  the  problem  from 
different  angles,  have  a  common  objective— the  analysis  of  aircraft  measured 
particle  data  to  yield  relations  between  liquid-water-content  and  some  other 
parameter  (radar  reflectivity,  temperature,  altitude)  which  may  then  be  applied 
to  determine  the  liquid  water  water  content  along  a  trajectory.  There  are  three 
possible  ways  of  determining  a  trajectory  M  in  the  absence  of  storm  specific  M-Z 
correlations.  In  the  first,  described  in  Section  2,  nominal  M-Z^  equations  were 
derived  from  one -second  aircraft  spectra  and  simultaneous  radar  reflectivity 
measurements  for  each  hydrometeor  type.  The  second  and  third  methods,  des¬ 
cribed  in  Section  3,  consisted  of  finding  relations  between  k  and  temperature  and 

0  5 

k  and  altitude,  and  using  the  relation  of  M  =  k  Z  . 

A  comparison  was  made  between  these  approaches  using  the  mean  values  of 
Zn,  temperature  and  altitude  of  the  hydrometeor  types  from  the  Kwajalein  data  in 

n 

Section  2.  The  mean  ZR  was  substituted  in  the  nominal  equation  to  give  one  esti¬ 
mate  of  M  for  each  type.  The  mean  temperature  was  used  with  the  Kwajalein 
specific  k  vs  temperature  relationship  to  give  a  second.  A  third  M  was  determined 
from  the  mean  altitude  and  the  k  vs  altitude  equation.  Table  5. 1  lists  the  percentage 
differences  that  result  in  the  derived  M  values  when  the  methods  from  Section  3 
are  compared  with  that  calculated  from  the  nominal  equations. 
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Table  5.  1  A  Comparison  of  M's  Derived  From 
the  k  Relationships  of  Temperature  and  Altitude 
at  Kwajalein  With  Those  Determined  From  the 
Nominal  Equations 


Hydrometeor 

Type 

Percentage 

AM/M  (Nominal  Eq) 
From 

k  vs  T  k  vs  Km 

Bullet  Rosettes 

6.  6 

7.7 

Small  Snow 

-25.3 

-23.6 

Large  Snow 

21.  0 

30.9 

Rain 

51.0 

30.9 

These  comparisons  highlight  the  variability  between  the  methods  using  pass 
averages  and  climatological  parameters  with  that  derived  from  the  one-second  air¬ 
craft  and  radar  measurements.  This  is  also  true  in  the  analysis  reported  in 
Section  4.  Although  the  data  used  in  Section  4  were  from  Wallops  Island,  and  less 
variable  than  those  from  the  tropical  convective  storms  of  Kwajalein,  the  limiting 
factor  in  defining  the  liquid  water  content  at  a  point  along  the  trajectory  was  still 
the  scatter  of  the  individual  estimates  around  the  pass  average.  The  results  of 
Section  4,  particularly  the  large  uncertainties  when  average  values  are  applied  in 
the  bullet  rosette  case,  demonstrate  that  situations  which  deviate  significantly  from 
climatological  averages  still  require  time -specific  M-Z  relations  derived  from 
simultaneous  aircraft-radar  measurements. 

Much  work  remains  to  be  done  ir.  this  field.  In  addition  to  an  adequate  theory 
of  how  the  melting  process  affects  radar  reflectivity,  mathematical  models  explain¬ 
ing  the  evolution  of  M-Z  relations  at  a  given  altitude  with  time  are  essential  before 
remote  measurements  of  M  will  be  possible.  As  data  become  available,  further 
studies  of  the  merits  of  using  M-Z  relations  tailored  for  particular  Z  profiles  will 
be  investigated. 
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LIST  OF  SYMBOLS  AND  ABBREVIATIONS 


a  power -function  coefficient 

b  power -function  exponent 

c  constant 

5C  degrees  Celsius 

d  minimum  drop  or  equivalent  melted  diameter 

D  drop  or  equivalent  melted  diameter 

D.  drop  or  equivalent  melted  diameter  of  classified  data 

maximum  drop  or  equivalent  melted  diameter 
D_j,  maximum  drop  or  equivalent  melted  diameter  that  is  able  to  be 

measured  with  a  particular  instrument 
H  altitude  in  kilometers 

i  index  of  summation 

k  aircraft  spectral  parameter 

kc  climatological  average  value  of  k 

L  measured  physical  size  of  ice  hydrometeors 

M  liquid  or  ice  water  content 

liquid  or  ice  water  content  from  aircraft  data 
Mc  liquid  or  ice  water  content  derived  from  cloud  probe 

Mc  climatological  liquid  of  ice  water  content 

Mp  distributed  liquid  or  ice  water  content 

estimated  liquid  or  ice  water  content  from  an  extrapolated  distribution 
Mp  liquid  or  ice  water  content  derived  from  precipitation  probe 


N 


D 

N. 


N„ 


M 


T 

Z 

ZA 

ZC 

Z 

c 

ZE 

Z 

P 

ZR 

AZ 

a 

A 


a> 


number  concentration  of  the  hydrometeor  drops  or  particles  per  m 
per  mm  bandwidth 
number  concentration  at  diameter  D 
number  concentration  for  classified  data 

the  "zero  intercept"  of  a  distribution  function  of  exponential  type 

3 

total  number  of  the  hydrometeor  particles  of  a  given  population  per  m 

number  of  terms  in  summation 

truncation  ratio  of  liquid  water  content 

truncation  ratio  of  number  concentration 

truncation  ratio  of  radar -reflectivity  factor 

temperature  in  3  Celsius 

radar-reflectivity  factor  for  Rayleigh  scattering 
radar-reflectivity  factor  computed  from  aircraft  data 
radar-reflectivity  factor  from  cloud  probe 
climatological  radar-reflectivity  factor 

estimated  radar-reflectivity  factor  from  an  extrapolated  distribution 
radar-reflectivity  factor  from  precipitation  probe 
radar-reflectivity  factor  as  measured  by  radar 

radar-reflectivity  factor  from  the  extrapolated  part  of  the  distribution 
coefficient  of  the  I  to  D  equation 
exponent  of  the  £  to  D  equation 

exponential  "slope  factor"  in  the  distribution  function  for  the  number 
concentration  of  the  drops  or  particles 
density  of  liquid  water  <g/cm  ) 

ratio  of  the  number  of  particles  in  class  i  to  the  total  number,  N^. 


AFCRL  Air  Force  Cambridge  Research  Laboratory,  Hanscom  AFB,  MA 

AFGL  Air  Force  Geophysics  Laboratory,  Hanscom  AFB,  MA 

BMO  Ballistic  Missile  Organization,  Norton  AFB,  CA 

KMR  Kwajalein  Missile  Range,  Kwajalein  Atoll,  Marshall  Islands 

NASA  National  Aeronautics  and  Space  Administration 

SAMSO  Space  and  Missile  Systems  Organization,  Los  Angeles  AFS,  CA 

WFTC  Wallops  Flight  Test  Center,  VA 

FMS  Particle  Measuring  Systems,  Boulder,  CO 

1- D  One -dimens ional  PMS  System 

2- D  Two-dimensional  PMS  System 

MIT  Massachusetts  Institute  of  Technology 

TRADEX  L  &  S  Band  Radar  used  at  KMR 
ALCOR  C -Band  Radar  used  at  KMR 
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